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Abstract: we derive exact analytical solutions of the Klein-Gordon equation for Makarov potential by means 
of the asymptotic iteration method. The energy eigenvalues are given in a closed form and the corresponding 
normalized eigenfunctions are obtained in terms of the generalized Laguerre polynomials and 
hypergeometrical functions . 

I. Introduction 

The Klein - Gordon equation (KG) is nowadays regarded as the relativistic form of the 
Schrodinger equation. It affords appropriate description for spin zero particles. Since the 
solution of the KG equation is often a complicated problem, use of mathematical methods is 
required. These include the variational method [1], the functional analysis method [2], the 
supersymetric approach [3], the Nikiforov - Uvarov method (NU) [4], the shifted 1/N 
expansion [5] and the asymptotic iteration method (AIM) [6]. The latter, however, is 
particularly suitable for obtaining solutions to such differential equations. Indeed, compared to 
other analytic techniques, the AIM offers an accurate and efficient investigation of the 
spectrum of many particles in relativistic and non relativistic quantum mechanics [7-10]. 

The KG equation has been investigated for several potentials as the generalized Woods- 
Saxon potential [11], Hulthen potential [12], pseudo and perturbed Coulomb potential [13-14] 
and equal scalar and vector potential [4]. In the present paper, we extend the investigation to 
Makarov potential [15] using the asymptotic iteration method [16]. The Makarov potential 
[15] offers an appropriate description of ring shaped molecules such as benzene CeH 6 and 
provides an interesting interaction model of deformed pairs of nuclei. 



II. Theory 

In spherical coordinates, the KG equation for a particle in general non-central 
potentialV(r, 0) can be written as 
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where E is the energy and M the mass of the particle. The non-central potential of Makarov 
[15] reads 

V (r i6 ) = -« + -£- + JpiL- (2) 
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where the first term represents the Coulomb potential while the second and the third are the 
short range ring shape terms. 

As usual we can write the wave function in the following form 

ip(r,9,<p) =^8(0)4>foO (3) 



then, the equation (1) will be separated into a radial part 

[^-^- + M) a - + E 2 - M 2 ]R(r) = (4) 
and angle dependent equations 
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where m and A, are separation constants. The solution of equation (5b) is well-known 

3>(<p) = Ae im<p ,(m = 0,±l,±2,...) (6) 

Only Eq. (4) and Eq.(5a) have to be solved. To this end, we use the asymptotic iteration 
method and the first step consists in the conversion of these equations to standard forms 
suitable to AIM applications [16]. So the radial equation (4) is rearranged as 

[iL_«£il> + 5£ + k2 ] fi(r) = ,7a) 

with 

X = t(t + l),s = ^j^- and k 2 = E 2 - M 2 (7b) 

and where the boundary condition is R(0) = 0. 

Now, the asymptotic behavior of the radial wave function suggests the following ansatz for 
R(r) 

R(r) = A{krY +1 e ikr f(r) (8) 

where A is a normalization constant. 

For this form of the wave function, the radial equation (7a) reads 

£f(r) = Ao(r) ^/(r) + s (r)/(r) (9) 

with 

Ao(r) = _p^) (l0a) 

5 (r) = -( 2ifer( 7 )+2s ) (10b) 
If, moreover, we introduce the new variablex = —2ikr, these equations transform into 



j^f(x) = Ao(x)j-f(x) + s (x)/(x) 



(11) 



l (x) = -(^) (12a) 

S (x) = (-^J (12b) 

According to the AIM procedure, the energy eigenvalues are then computed by means of the 
following termination condition [16] 

8 = s n A n _ t - X n s n _ x = (13) 
for a given x > 0, with the sequences 

= CiW + s„_i(x) + A W4-iW (14a) 

s n 00 = Sn-iCO + s (x)l n _ 1 (x), n= 1, 2, 3, ... (14b) 
For n = 10 iterations, the obtained solutions are 



In general form, we have 



k = -^-,k = -^-,k = -^-,k = -^-, (15) 
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k = -^—,N=0,l,2, (16) 
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Substituting k and s by their expressions given in Eq. (7b), we finally derive the exact 
eingenvalues E N of the radial part of the Klein Gordon equation 

(£+l+N) 2 -4- 

E N =M- — £ (17) 
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where N is the radial quantum number. 

Now we turn solutions of the polar angle equation (5 a). By introducing a new variable 

y = - (1 + cos9), it becomes 

[y(i - y) |-y(i -y)± + W + i)ytt - y) - a 2 - Q> 2 - a 2 ) y ] c(y) = o (is) 

with 

a = Jm 2 +p' 2 - Y ' 2 b _ Jm 2 +P' 2 + Y " z ^ 

fi' = (E + M)/? and y = (E + M)y. Since G(y) must satisfy the boundary conditions 
G(0) = G(l) = 0, we use the following ansatz 



G(y) = By a (l-y) b f{y) (20) 
Substitution of this wave function into Eq. (18) leads to the AIM standard forms 

£sf(y) = AoGO |r/(y) + s (y)f(y) (21) 
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A oGO = (22a) 



s o (y) = — } (22b) 

We may then derive the eigenvalues of Eq. (18) by considering AIM equation (12) for a given 
y > 1. For 10 iterations, we obtain 

£ = a + b,£ = a + b + 1, £ = a + b + 2, £ = a + b + 3, ... (22) 

Or equivatently 

£ = a + b + n, n = 0, 1, 2, 3, (23) 

Substituting a and b by their expressions given in (19) leads to 



t n = f 2+ e'w + l£±l + n (24) 

Once this condition for £ is inserted in the previously derived expression (17) we get the exact 
energy eigenvalues for a bound particle in the Makarov potential 
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The obtained energy spectrum is identical to the one determined in the paper [4] via the 
Nikiforov - Uvanov method. 

Next, we proceed to look for the radial wave function R(r). For this, we have to determine 
the function f(x) by solving the differential equation (11). Substituting k by its expression given 
in Eq. (16), this equation transforms into Kummer equation 

^ /W + (^)A /W + W/(X)=0 (26) 

whose solutions as x — > are the confluent hypergeometrical functions, 

fix) = C N F(-N, 2£ + 2, x) (27) 



Thanks to the relation between Kummer functions F (a, b; x) and the generalized Laguerre 
polynomials L^f +1 (x), the radial wave function can be written 

R«W = D N£ x <+1 e-h 2 * +1 (x) (28) 

where D N{ is a normalization constant computed via the orthogonality relation of Laguerre 
polynomials [17] 

, _ r (£+M) g (n'-l-l)! 1 1/2 nQ . 
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where n' = iV + ^ + 1 and finally we derive the expression of the radial wave function 

r (g+MMn'-l-l)! l 1/2 / (g+M)« ^ +1 , +1 JgffiHr 2 < +1 r (E+M) g 

Similarly, the angular wave function G(y) is obtained by solving the hypergeometric 
differential equation whose solutions are given by [17] 

/GO = 2 Fd-n,a + b + £ + l;2a + l;y) (31) 

where 2F1 is a special case of the generalized hypergeometric function. Therefore, by 
combining Eqs. (20) and (31), the polar angle wave functions are found to be 

Gemiy) = iW a U - jf 2 Fd-n, a + b + £ + 1; 2a + l;y) (32) 

From the orthonormality of hypergeometric functions, we determined the normalization 
constant 

„ 1 J T(l+q+fc+l)r(l+a-fc+l)(2l+l) j 1 / 2 

{m ~ T(2a+1) L 2(t-a-b)\Y(t-a+b+l) \ ^ ' 

where the parameters a and b are provided by Eq.(19). 
III. Conclusion 

In the present paper, the bound states of the Klein Gordon equation for Makarov potential 
are investigated. Using the asymptotic iteration method, we have derived the exact expression 
of the energy eigenvalues and the corresponding normalized eigenfunctions in terms of the 
generalized Laguerre polynomials and hypergeometrical functions. The obtained analytical 
results are very useful in nuclear physics and quantum chemistry. 

Acknowledgments: 

This study is partially supported by the CNRST Moroccan research program PROTAS III, 
D16/04. 



IV. References 



[I] B. Batiha, Australian Journal of Basic and Applied Sciences, 3(4): 3876-3890 
( 2009) 

[2] Ying Zhang, Phys. Scr. 78 015006 (2008) 

[3] Gang Chen, Zi-Dong Chen and Zhi-Mei Lou, Phys. Letters A 
Volume 331, Issue 6, (2004)374-377 

[4] Fevziye Yasuk, Aysen Durmus and Ismail Boztosun, J. Math. Phys. 47, 082302 
(2006) 

[5] Mustafa, Omar; Sever, Ramazan, Phys. Rev. A, vol. 44, no7, pp. 4142-4144(1991) 
[6] Eser Olgar 2008 Chinese Phys. Lett. 25 1939 

[7] Brodie Champion, Richard L. Hall, Nasser Saad, arXiv:0802.2072vl [math-ph] 
[8] Inci, I. ; Boztosun, I. and Bonatsos, D. AIP Conference Proceedings V.1072 (2008) 
[9] Soylu and I. Boztosun, Physica E Volume 40, Issue 3, January 2008, pp. 443-448 
[10] I. Boztosun, M. Karakoc, F. Yasuk and A. Durmus, arXiv:math-ph/0604040vl 

[II] Ikhdair, S. M. Sever, R, ANNALEN DER PHYSIK -LEIPZIG-, VOL 16; 
NUMBER 3, pages 218-232 (2007) 

[12] SIMSEKMehmet and FGRIFES Harun, J. Phys. A, vol. 37, n°15, pp. 4379-4393 
(2004) 

[13] Xiao Yong-Jun and Long Zheng- Wen, Commun. Theor. Phys. 53 54 (2010) 

[14] T. Barakat, Annals of Physics, V. 324, Issue 3,2009, pp. 725-733 

[15] Makarov A. A., et al., 1967, Nuovo Cimento A 52 1061. 

[16] H. Ciftci, RL. Hall andN. Saad, J.Phys.A 36(2003)11807-11816 

[17] Chang- Yuan Chen, Cheng-Lin Liu and Fa-Lin Lu, Phys. Lett. A 374(2010)1346- 

1349. 



